TUW-04-19 



Collisional Energy Loss of a Heavy Quark in an Anisotropic 

Quark-Gluon Plasma 

Paul Romatschke and Michael Strickland 

Institut fiir Theoretische Physik, Technische Universitdt Wien, 
Wiedner Hauptstrasse 8-10, A-IO4O Vienna, Austria 



Abstract 

We compute the leading-order cohisional energy loss of a heavy quark propagating through 
a quark-gluon plasma in which the quark and gluon distributions are anisotropic in momentum 
space. Following the calculation outlined for QED in an earlier work we indicate the differences 
encountered in QCD and their effect on the collisional energy loss results. For a 20 GeV bottom 
quark we show that momentum space anisotropics can result in the collisional heavy quark energy 
loss varying with the angle of propagation by up to 50%. For low velocity quarks we show that 
anisotropics result in energy gain instead of energy loss with the energy gain focused in such a way 
as to accelerate particles along the anisotropy direction thereby reducing the momentum-space 
anisotropy. The origin of this negative energy loss is explicitly identified as being related to the 
presence of plasma instabilities in the system. 

PACS numbers: 11.15.Bt, 04.25.Nx, ll.lO.Wx, 12.38.Mh 
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I. INTRODUCTION 



An understanding of the production, propagation, and hadronization of heavy quarks in rel- 
ativistic heavy-ion cohisions is important for predicting a number of experimental observables 
including the heavy-meson spectrum, the single lepton spectrum, and the dilepton spectrum. The 
first experimental results for the inclusive electron spectrum have been reported ^ in addition to 
the first measurements of J/ ip production at RHIC j^l ■ The measurement of the inclusive electron 
spectrum allows for a determination of heavy quark energy loss since it is primarily due to the 
semi-leptonic decay of charm quarks. The heavy quark energy loss comes into play since it is 
necessary in order to predict the heavy quark energy at the decay point. It is therefore important 
to have a thorough theoretical understanding of heavy quark energy loss for a proper comparison 
with the experimental results. In this paper we will show that in QCD there is a modification of 
the leading-order (collisional) heavy quark energy loss if there is a momentum-space anisotropy in 
the underlying quark and gluon distribution functions. 

Here we will focus on the collisional energy loss using the techniques of Braaten and Thoma 
Q which have been extended to anisotropic systems in Ref. [H]- For a review of the radiative 
energy loss in isotropic systems see Refs. p and The method employed by Braaten and 
Thoma gives the complete result for the collisional energy loss by considering independently the 
contributions from soft (involving momenta (7 ~ m^)) and hard {q ~ T) gluon exchange. The two 
scales, mo and T, are then separated by an arbitrary momentum scale q* which cuts off the UV 
and IR divergences in the soft and hard contributions, respectively. Moreover, it was found that 
in the weak-coupling or high-temperature limit <C 1 the condition mo ^ q* ^ T could be used 
to expand the resulting integral expressions for the soft and hard contributions further, giving an 
analytic result for the energy loss independent of the separation scale q* However, in the 

case of QCD at temperatures accessible by heavy-ion collision experiments the coupling constant 
is quite large (a^ ~ 0.2 — 0.3) and as a consequence the analytic expression of Braaten and Thoma 
can give unphysical results for the energy loss as we will demonstrate. 

On the good news side, it turns out that when one does not expand the integral expressions 
with respect to the condition m/) <C <C T the isotropic contributions to the energy loss always 
stay positive, even for very large coupling, as has been shown for QED |^. However, this comes 
at the expense of giving up independence of the complete result on the scale q*, which then has to 
be fixed somehow. Fortunately, it turns out that in the weak-coupling limit this scale dependence 
becomes very small and as a consequence, unless q* is taken to be very large {q* ^> T) or very 
small {q* <C ttid), one recovers the analytic results found by Braaten and Thoma. When the 
coupling is increased one can then still fix q* by the "principle of minimal sensitivity", which 
means that q* is chosen such that the energy loss (and therefore also the variation with respect to 
q*) is minimized. In addition, the principle of minimal sensitivity provides a way to estimate the 
theoretical uncertainty of the result by varying q* by a fixed amount around the point at which 
the energy loss is least sensitive to this scale. 

The main reason why a calculation of this t ype is important to perform is that the presence of 
plasma instabilities 0, 0, 0, 0, 11 , 0, 0, 0? 0? Q| iii anisotropic systems could have 

a significant effect on observables like the heavy quark energy loss. In fact, the presence of such 
instabilities naively renders the calculation of the soft part divergent; however, there is protection 
mechanism dubbed "dynamical shielding" which renders the collisional energy loss finite for QED 
This protection mechanism trivially extends to the case of QCD since the only thing that 
changes in going from soft QED to soft QCD at leading-order is the numerical value of the Debye 
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mass. However, at large coupling the presence of instabilities and associated poles on unphysical 
Riemann sheets (l^ causes significant changes in the soft energy loss contribution for both QED 
and QCD. In fact, as we will discuss, the unphysical poles can even change the sign of the heavy 
quark energy loss at low momentum turning it instead into energy gain. 

The paper is organized as follows: In Sec. we review the calculation of the collisional energy 
loss including the case of very strong anisotropics. In Sec. lIIII we discuss the role that modes on the 
unphysical sheets play in the calculation. In Sec. IIVI we present results for the collisional energy 
loss of a heavy quark for isotropic and anisotropic systems. We give our conclusions in Sec. IVl 



II. COLLISIONAL ENERGY LOSS IN QCD 

Following Refs. [^jSIHjHI the collisional energy loss of a heavy quark with velocity v is divided 
into a soft contribution and a hard contribution. The soft contribution within QCD has exactly 
the same form as in QED except that the expression for the Debye mass is different and there is an 
overall Casimir which has to be taken into account. Therefore the soft contribution from Ref. |^ 
can be taken directly and only needs to be rescaled appropriately. The hard contribution in QCD is, 
however, different for two reasons: a new gluon scattering diagram appears due to the three-gluon 
coupling and the graphs which are topologically equivalent to the QED fermion-photon graphs no 
longer cancel due to non-trivial Casimir invariants. In the next two subsections we present explicit 
integral expressions for both the soft and hard contributions to the energy loss in the limit of 
infinite quark mass. In the final subsection we discuss how the limit of infinite anisotropy can be 
taken. 



A. Soft contribution 

The soft contribution in the case of QCD is given by the expression 

- = 2^ '-/Iff <- ^ [^"«) - ^»'<«1 • 

where Q is the quark color charge in the fundamental representation with = g'^{N^ — l)/{2Nc) 
and the two propagators A*-'((5) and Aq (Q) denote the hard-loop and free gluon propagators, 
respectively. The hard-loop gluon propagator in an anisotropic system has been calculated in 
Ref. and can be expressed as 

A{Q) = Aa{Q) [A - C] + Ag{Q) [{q^ -u;^ + a + j)B + {p- u;^)C - 5D] , (2) 

where the tensor basis for the spacelike components of a system with one preferred direction 
specified by n is: A^^ = 5^^ - W /k"^, B'^ = k'^/k"^, C^^ = n'n^/h^, D'^ = k'n^ + k^n' with 
= A^^ . The propagators A a and Aq are then given by 

^aHQ) = q^-oJ^ + a, 

^gHQ) = (g' - + a + 7)(/3 - u;^) - q^n^S^ , (3) 



with a, P, 7, and 5 being the hard-loop gluon structure functions. Except for some special cases 
(certain angles of propagation, infinite anisotropy, etc.), the analytic expressions for the structure 
functions are quite complicated in form and it is more convenient to use a numerical evaluation of 
their integral representations 
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The structure functions of anisotropic systems depend on the precise form of the quark and gluon 
distribution functions. Here we will assume that the distribution function is spatially homogeneous 
and therefore given by 

/(p) ^ 2Nf (n(p) + n(p)) + AN,ng{p) , (4) 

where n, n, and Ug are the distribution functions of quarks, anti-quarks, and gluons, respectively.^ 
We will further assume that the anisotropic distribution function, /(p), can be obtained from an 
isotropic distribution function by rescaling one direction in momentum space 



/(P) = /iso + ^(P • n)' 



(5) 



where N{^) = -v/T+T is a normalization constant with ^ the strength of the anisotropy. The 
anisotropy strength ^ is in the range — 1 < ^ < oo with ^ = corresponding to the original 
isotropic distribution. In order to simplify the analysis we take quarks and gluons to have the 
same anisotropy parameter although in general they can be different. Even though its effects 
might be interesting, different anisotropics for quarks and gluons will not be studied in this paper. 

The soft contribution to the collisional energy loss of a heavy quark in an anisotropic quark- 
gluon plasma is then given by 
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(8) 



and a) = q • V, and 5 = q6. In order to regulate the soft contribution it is necessary to introduce a 
UV momentum cutoff q* on the q integration. 

Note that in the case of an anisotropic system unstable modes are present that manifest them- 
selves as potentially unregulated poles of the propagator in the static limit. However, as discussed 
in our previous paper Q the mechanism of "dynamical shielding" protects the collisional energy 



^ There is a factor two difference between Eq. (0} and the equivalent expression given in Ref. 5]. We have 
changed our notation in this work so that all symmetry factors are specified in /(p). Final results are 
unaffected by this redefinition. 
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loss from these potentially devastating singularities in the soft contribution. These singularities 
could arise in terms containing, for example, {q^ + due to the fact that in the static limit 

the structure function a is negative- valued. However, if one takes the static limit of a carefully we 
find 

\{nia{uj,q) = M^{-l + iDu) + 0{uj'^) , (9) 

where M and D depend on the angle of propagation with respect to the anisotropy vector and the 
strength of the anisotropy. As long as D is non-vanishing the singularities are regularized because 
of the u) in the numerator of Eq. and we call the singularity "dynamically shielded" . A proof 
of dynamical shielding for very weak and strong anisotropics can be found in Refs. [H] an d ll^ . 

In the isotropic limit = 0) Eq. Q corresponds to the result of Thoma and Gyulassy 19] and 
one obtains the soft contribution found by Braaten and Thoma 0] when expanding Eq. (jH)) under 
the additional assumption q* ^ niD • The effect of this expansion is that the result of Ref . 0| 
becomes negative for small q* /niD (but allows the calculation of an analytic result for the isotropic 
collisional energy loss) while the unexpanded result is positive for all q* /uid- 



B. Hard Contribution 



The hard contribution can be separated into two parts: one contribution coming from the 
scattering of the heavy quark on quarks in the plasma and another that takes into account the 
scattering on plasma gluons (corresponding to the tree- level diagrams shown in Fig. ^) . Assuming 
the velocity of the quark to be much higher then the ratio of the plasma temperature to the energy 
of the quark, v ^ T/E, the contribution coming from quark-quark scattering can be reduced to 
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after performing the Dirac traces and evaluating the sum over spins. The contribution coming 
from quark-gluon scattering gives 

_ (4^)3q2 f d?k ng{k) f d^k' l + ng{k') 
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Here n(k) and ng(k) are the anisotropic versions of the tree-level Fermi-Dirac and Bose-Einstein 
distribution functions at zero chemical potential specified in Eq. Q), to = k' — k, and q = k' — k. 
Note also that q* acts as an IR cutoff for the q integration. Since the integrand is odd under the 
interchange k k' the terms involving the products re(k)n(k') and ng(k)ng(k') vanish since they 
are symmetric. Following Ref. 5], some of the integrations can be performed giving 

_/clW^\^« _8alNf{q*fT^^/TT^ ''^ ^ .2^ 
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FIG. 1: Tree- level Feynman diagrams for the scattering processes Qq — > Qq (first diagram) and 
Qg ^ Qg (remaining diagrams). 
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where T denotes the unwieldy expression 

G(42;2 sin^ 9^ siiP 9g - {1 - cos^ 6'<y + 2 cos 6'gz(cos 9^ - w))^) 
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where x = g*2;\/l + ^(na; sin6'fe cos 0^ + cos 9kY /T, with = n • v and 1 = + n^. The 
remaining integrations have to be performed numerically. 

Similar to what has been found for the soft part, in the isotropic limit the hard contribution 
becomes that of Braaten and Thoma |J] when expanding Eq. ()12() with respect to q* <^ T; also, the 
expanded result becomes negative for large q* /T, while the unexpanded result is always positive 
in the isotropic limit. 



C. Very strong anisotropies 

There are however a few special cases where the structure functions become quite simple, notably 
the extreme anisotropy limit, oo, which corresponds to the parton distribution function having 
the form 

/(p) = F{pMPz) ■ (15) 

The resulting expressions for the structure functions have been calculated in Ref. and this 
special case has been treated in more detail in a separate publication 0|. Using these structure 
functions one is able to calculate the soft contribution to the energy loss in the large ^ case by 
evaluation of Eq. However, in the extreme limit — > oo this is complicated by the presence of 
troublesome spacelike quasiparticle modes For large but finite ^ these modes are integrable 
because they move off the physical sheet onto the neighboring unphysical sheets but for ^ = oo they 
result in a pole in the integrand that needs to be regulated. Fortunately, a natural regularization 
is provided by the fact that for any finite ^, these modes are outside of the physical spectrum; 
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therefore, when keeping ^ finite and sending it to infinity only after having done the integrations, 
one is able to produce well-defined results.^ 

The hard contribution in the large ^ case is easily found by using Eq. H12|) together with 

hm y/l + (,Fi {q*z\l l + i{k- hf^ = 5{k ■ h) J dy Fi{q* I + y"^) (16) 

and similarly for F2- Using 

k ■ n = Ux sin 9^ cos + Uz cos 9k , (17) 
the delta function can be cast into the form 

6{k . n) = ^ ^^'^^ - '^"^ , (18) 



with 



ni sin 9k — nt cos^ 9k 



Ux cos </>o = —nzCoi9k ■ (19) 



Using the symmetry In — (p ^ (p and substituting cos 9k — > Hx cos 9k the remaining integrations 
are easily implemented. 

III. EFFECTS OF MODES ON THE UNPHYSICAL SHEET 

In a previous paper we have discussed the presence of nearly-spacelike "unphysical modes" in 
anisotropic plasmas ^3]- In that work we showed that for large anisotropics new relevant singu- 
larities are present on neighboring unphysical Riemann sheets. The presence of these unphysical 
poles has a resonance-like effect on propagators on the physical sheet and therefore can influence 
observables which are sensitive to spacelike momentum. Before presenting our numerical results 
we would like to first discuss how the presence of these unphysical modes affects the soft contri- 
bution to the energy loss. In order to do this as clearly as possible we consider the case when the 
heavy quark is propagating along the anisotropy direction, i.e. 6^ = 0, in which case the azimuthal 
integration in Eq. becomes trivial since a) — > v cos 9q. In this case the soft contribution to the 
energy loss takes the form 

fdW\ T r\ Q^sm9g q^A + B 

- I J so...^0 = Wrv'^Jo Jo q^C + q-V + 8 ' 

with "n? = n • V = 1 — cos^ 6q understood in Eq. (jSJ. In this special case, we can easily relate 
Eq. H^Uj) to the collective modes in the system: this is done by investigating the integrand at fixed 
angles 9q, which represents the contribution from the straight line oj = qvcos9q from g = to 
q = q* ■ This integration path can then be directly compared to the unphysical mode dispersion 
relation cuj = uji{q,9q) where i = {A,B} indexes the unphysical mode in question. 

As a further simplification we next consider the limit ^ ^ oo for which the dispersion relations 
for the collective modes on the unphysical sheets are known for arbitrary angles 9q of propagation 
[l^ . This limit is convenient because it is possible to analytically determine the structure functions 
in this case, however, the general conclusions discussed below still apply when ^ is finite. In the 
top plots of Fig. I^Jwe show the integration paths corresponding to (a) v = 0.5 and (b) v = 0.99 as 



^ In practice it is more convenient to use the ^ — > cx) expressions for the structure functions and regulate 
the integrand "by hand" so that it agrees with the finite £, case. 
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q/niD 



FIG. 2: In the top plots we show the integration paths corresponding to (a) v = 0.5 and (b) 
V = 0.99 as solid lines along with the corresponding dispersion relations for the two infinite-^ 
unphysical modes (dashed and dotted lines labelled as "Mode A" and "Mode B"). In the bottom 
plots we show the corresponding value of the integrand of Eq. (|20j) evaluated along the integration 
paths shown in the top plots. We also plot the integrand of Eq. (PH]) for ^ = {1000,100,10} as 
decreasingly lighter gray lines in the bottom plot. Note that in both (a) and (b) the integration 
paths in the top plots should be understood to terminate at q = q* . 



solid lines along with the corresponding dispersion relations for the two infinite-^ unphysical modes 
(dashed and dotted lines labelled as "Mode A" and "Mode B"). In the bottom plots of Fig. [21 we 
show the corresponding value of the integrand of Eq. ()20() evaluated along the integration paths 
shown in the top plots. Note that in both (a) and (b) the integration paths should be understood 
to terminate &t q = q* . 

Let us discuss the case oiv = 0.5 first: as indicated by the vertical line in Fig.[2fa), the extremal 
contribution to the energy loss occurs when the integration path comes closest to the mode A. 
However, it turns out the contribution from this mode is negative, representing an inversion of 
the usual Landau-damping mechanism (we will discuss this issue later on). This situation is not 
limited to infinite anisotropies. In the bottom plot of Fig. I2ta) we show the integrand of Eq. H20|l 
for ^ = {1000,100,10} as decreasingly lighter gray lines. As we can see the effect remains but 
becomes less as the anisotropy parameter is decreased. 

In Fig. Hfb) we show the same quantities but for v = 0.99.^ As indicated by the vertical line 
in Fig. EKb), the extremal contribution to the energy loss in this case occurs when the integration 



^ There is no "regular" contribution to the soft e nergy l oss for ^ — > oo, since the structure functions have 
non- vanishing imaginary part only for uj < sin 6g |l6lll7l| . which together with uj — v cos 9q gives v < tan 9q. 
However, since mode B moves onto the physical sheet for ^ ^ oo, the argument of the logarithm in Eq. ©, 
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path comes closest to the mode B. However, in contrast, mode B contributions to the integrand 
are positive. Again, this behavior persists at finite as shown by the corresponding plots of the 
integrand of Eq. ()2U|) for ^ = {1000, 100, 10} in Fig.EJb) (shown as increasingly lighter gray lines). 
However, the peak of the integrand for mode B contributions moves considerably more than the 
mode A contributions as ^ is decreased. 

As shown by the light gray lines in the lower plots of Fig. |21 the unphysical A and B modes 
seem to have the same qualitative effect on the energy loss at finite ^ as they do in the limit of 
^ ^ oo. However, since the modes move farther away from the physical region as ^ is decreased, 
it is clear that their influence becomes smaller, as already argued in Ref. [13]- Nevertheless, it is 
interesting to note that even for small values of ^, mode A is capable of driving the soft energy loss 
contribution negative for small velocities v. Indeed, it turns out that in the limit v ^ 0, one finds 
that the soft energy loss is purely negative (independent of the angle), as long as ^ > 0. 

This negative energy loss contribution is, however, not due to the fact that a particle with 
sub-thermal velocities should receive energy from the medium, since our approximation was based 
on the assumption of an infinitely heavy quark (see the discussion about limitations in the next 
section). Rather, since mode A is the analytic continuation of the plasma instability present on 
the physical sheet, we interpret the gain of energy as being intimately related to the presence of 
such instabilities. This is not a surprising finding given that if there are truly unstable modes in 
the system the energy for their growth must come from a mechanism which is also collective since 
it would be hard to imagine that a collisional mechanism could ever be fast enough. 

The mode A contributions to the soft energy "loss" represent such a collective mechanism for 
the efficient transfer of energy from hard to soft scales. As we will show below the energy gain 
is focused in such a way that particles with small velocities are accelerated out-of-plane thereby 
reducing the anisotropy. In closing we note that for angles different than 9^ = 0, the connection 
between the modes on the neighboring unphysical sheets and the integrand of the energy loss 
remains but it is more difficult to disentangle since there is also an averaging over the azimuthal 
angle as well. 



IV. HEAVY QUARK ENERGY LOSS 

The collisional energy loss of a heavy quark in an anisotropic quark-gluon plasma to leading 
order in the coupling is given by adding the contributions Eq. @ and Eq. ((T^ 
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In all results presented below we will assume A''^ = 3 and Nf = 2. 



A. Isotropic results 

In order to compare with existing calculations we first consider the isotropic limit (^ = 0). In 
this case the collisional energy loss is a function of the strong coupling as, the particle velocity v, 
the temperature T, and the momentum separation scale q*. The q* dependence of the result is 



though real, changes sign there, picking up a phase which results in the (5-function like peak in the lower 
plot of Fig.|a;b) (black line). 
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FIG. 3: Isotropic energy loss of a bottom (a) and charm quark (b) as a function of momentum p 
for Us = 0.3. Shown are the respective results from Bjorken |20| (dotted line) and Braaten and 
Thoma 0] (dashed line) as well as the evaluation of Eq. (PT|) (full line) with a variation of q* (gray 
band). 



found to become weak for small values of the coupling limit (similar to what has been found in 
QED 5]), corresponding to the original result by Braaten and Thoma 4]; for larger values of the 
coupling one fixes q* = qP™^ using the principle of minimum sensitivity 

d fdW\ 
dq* \ dx J „*=„• 



. (22) 



Note that in this case — {dW/dx)\g*^g, always serves as a lower bound on the result for the 
energy loss. To get an estimate of how strongly the result depends on this special value of q*, one 
can e.g. vary q^^^ by a certain factor Cg* and evaluate the energy loss at QpmsCg* and Qpms/cg*- 
This factor Cg* is in principle arbitrary; here we fix it to be Cg* = 2 to be consistent with what has 
been done in QED [a]. 

In Fig. 121 various results for the heavy quark energy loss in an isotropic quark-gluon plasma are 
compared: shown are the results from Bjorken Braaten and Thoma ^, as well as Eq. (|21jl 
with ^ = at g* = q^^^g together with its variation using Cg* = 2. Both plots assume = 0.3. 
Fig. Ol^a) shows the energy loss of a bottom quark with mass Mq = 5 GeV while in Fig. Ofb) the 
energy loss of a charm quark with mass Mq = 1.5 GeV is plotted, both as a function of their 
momenta 

In closing we would like to emphasize that our result is not only compatible with previous 
results on isotropic collisional energy loss but, in fact, represents the most complete calculation of 
this quantity to date. In contrast to the result of Bjorken it gives the correct leading-order 
result in the weak-coupling limit as does the result of Braaten and Thoma 0]; however, as noted 
before, at realistic values of the coupling and small momentum the Braaten and Thoma result 
breaks down completely. In the isotropic case our result is positive definite as it should be for an 
infinitely massive quark. In addition the use of a principle-of-minimal sensitivity to fix the scale 
gives the added benefit of being able to place theoretical error bars on the results obtained. 
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FIG. 4: Sketch of the vectors relevant to the dh'ectional dependence of the energy loss: n is the 
direction of the anisotropy, v the velocity of the quark, and the angle between these. We have 
also plotted an isocontour of a > distribution function for comparison. 



B. Limitations 

Since Eq. H21|) has been derived assuming an infinitely heavy quark we have to determine what 
will happen if we instead have only a heavy quark like a cha rm or bottom quark. The chief 
consequence is that for sub-thermal velocities defined by f < ^Jt/Mq the approximation breaks 
down because quarks with those velocities will gain energy from collisions with other particles in 
the heat bath instead of losing energy. Below this threshold the energy loss would then turn into 
energy gain. A semi-quantitative estimate of the velocity where the isotropic energy loss becomes 
negative has been performed in Ref. by repeating the above calculation in the limit v — > and 



for weak couplings, finding v ^ ^3T/Mq. Using T = 250 MeV and the masses above, this velocity 
corresponds to p = 2.1 GeV and p = 1.5 GeV for bottom and charm, respectively. Similarly, 
Eq. (|21|) also breaks down for ultrarelativistic energies E > M^/T, with the cross-over energy 
having been determined to be -Ecross — 1-8Mq/T (corresponding to v > 0.995 for both charm and 
bottom quarks). 

In contrast the Braaten and Thoma results shown in Fig. |H1 turn negative at a momentum of 
p ~ 5.7 GeV for the bottom and p ~ 1.71 GeV for the charm quark. This behavior is not due 
to the physical reason discussed above but rather due to a failure of the extrapolation from the 
weak-coupling limit to realistic couplings For the unexpanded isotropic result Eq. (|2H) . this 
unphysical behavior does not occur (it is always positive) and one can therefore expect it to be 



valid for velocities down to the original estimate v ~ yST/Mg. In anisotropic systems, however, 
this rule of thumb no longer applies and negative energy loss does not require the heavy quark 
to have sub-thermal velocity. Instead, as already noted by Lifshitz and Pitaevskii this type 
of energy gain is connected to the instabilities of the system, as we have argued in the previous 
section. 

Before moving on to the anisotropic results let us briefly comment on the variation estimate 
based on setting Cq* =2: while for QED this represents a very conservative choice on the variations 
since m/j <C qpms ^ 27rT, in QCD one encounters situations with Qp^s > ^vrT or g*,^^ < m/) for 
realistic couplings like = 0.3 since niD and 2ttT are no longer well separated. Strictly speaking, 
the method we are using is breaking down at realistic couplings; nevertheless, by using variations 
with Cq* = 2 we can hope to obtain reasonable estimates of the leading-order results when the 
difference between these variations is not too large. Finally, it should be noted that no estimate 
of the next-to-leading order (NLO) radiative corrections to the energy loss has been made here. 
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FIG. 5: Anisotropic energy loss for ^ = 10.0, 
= 1 X 10^^, V = 0.1, and 6*^, = as a 
function of q* . Dotted, dashed, and solid 
lines are the hard contribution, the soft con- 
tribution, and their sum, respectively. 
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FIG. 6: Anisotropic energy loss for ^ = 10.0, 
a<j = 1 X 10^^, V = 0.1, and 0^ = as a 
function of q* . Dotted, dashed, and solid 
lines are the hard contribution, the soft con- 
tribution, and their sum, respectively. 



Therefore, it should be kept in mind that for QCD with large realistic coupling the inclusion of 
these NLO corrections might give energy loss results that are not covered by the variations of q* 
in the leading-order result, so these variations should be interpreted with care. 



C. Anisotropic results 

For an anisotropic system with (gluon and quark) anisotropy strength ^, the collisional energy 
loss result Eq. (PT|) depends on q* , T, and v as was the case for the isotropic energy loss, but in gen- 
eral also on the angle of the particle propagation with respect to the anisotropy vector cos 0„ = v • n 
(sketched in Fig. |1J) . As in the isotropic case the total result is obtained by combining the soft and 
hard contributions to the energy loss. In Figs. I5I8I we plot the hard and soft contributions to the 
energy loss together with their sum as a function of q* /T for = {1 x 10~^, 1 x 10~^, 0.01, 0.3} 
corresponding to m^jT = {0.041,0.13,0.41,2.2}. As can be seen from these Figures the soft con- 
tribution (dashed line) is always negative over some range of q* /T. For small couplings (see Fig. [5J 
this negative contribution is in a region where there is always a larger positive hard contribution. 
However, as the coupling constant is increased we see that the range of q* /T over which the soft 
contribution is negative moves to higher scales and eventually dominates the sum of the soft and 
hard contributions over a range of q* (see Fig. |SJ . 

For all ^ we find that the soft energy loss has a minimum at some value of q* . For small 
V this minimum can even be negative as explicitly demonstrated in Figs. [7| and |S1 Considering 
again 9^ = {) for guidance we notice that the limits on the integration in Eq. ()20() correspond to a 
triangular region defined hy Q < to < qv and q < q*. From this we observe that when q* is small 
we only receive contributions from mode A which are predominantly negative. As q* is increased, 
however, we start to receive large positive contributions from mode B. Defining q* as the point 
where the upper boundary defined hy oj = qv hits mode B we obtain q* = ujb/v. In order to 
estimate the value of q* we again consider the limit £ — > oo. The dispersion relation for mode B 



at 9q = and £ = cxd is known and is g-independent 17 1 



= ■ (24) 
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FIG. 7: Anisotropic energy loss for ^ = 10.0, 
Us = 0.01, V = 0.1, and 0^, = as a function 
of q* . Dotted, dashed, and solid lines are the 
hard contribution, the soft contribution, and 
their sum, respectively. 
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FIG. 8: Anisotropic energy loss for ^ = 10.0, 
a<j = 0.3, V = 0.1, and ^t, = as a function 
of q* . Dotted, dashed, and solid lines are the 
hard contribution, the soft contribution, and 
their sum, respectively. 



Using this we obtain 

For larger values of q* , the negative contributions also become larger, but there is also now a 
strongly positive contribution from mode B (even for nonzero 6q), so that the overall result turns 
out to be more positive than at q*. On the other hand, for q* < q*, the main negative contributions 
become smaller since the "region" they are integrated over has shrunk. Therefore, the minimum 
of the soft contribution to the energy loss for ^ ^ oo at 0^, = is obtained around q*. This can be 
used as a rough guideline for the location of this minimum also for large but finite ^, however, as 
^ is decreased q* moves to lower values of q* and the depth of the minimum decreases. 

Returning to the results for the total anisotropic energy loss (hard plus soft) we again fix the 
scale q* by the application of the principle of minimal sensitivity 1)22(1 as in the isotropic case; 
however, in order to keep the figures presented as understandable as possible we do not show the 
corresponding variation of q* by Cq* =2. Instead we note that the sensitivity of the anisotropic 
results to this scale is similar to the isotropic variation indicated by the gray band in Fig. ^ 
however, there is an increase in the magnitude of the variation as ^ is increased. For example, 
at ^ = oo the corresponding variation is approximately three times the corresponding isotropic 
variation. 

In Fig. El we show the energy loss Eq. ()21() for a heavy quark propagating parallel to the 
anisotropy direction (6^ = 0) and perpendicular to the anisotropy direction (By = tt/2) as a 
function of the heavy quark velocity v with = 0.3 and ^ = 1. The results shown are normalized 
to the isotropic energy loss. As can be seen from this Figure for small velocities the energy loss is 
larger at 9^ = Tr/2 while for large velocities it is larger at 9^ = 0. This is similar to the behavior 
found for a QED plasma 

In Fig. ^1 we show the corresponding energy loss for a bottom quark with Mq = 5 GeV using 
the relation between velocity and momentum given in Eq. (|23|1 . We have not normalized this result 
as before but instead include the corresponding isotropic result as a solid line for comparison. As 
can be seen from this Figure for ^ = 1 and = 0.3 the energy loss of a bottom quark with 

= vr/2 is almost indistinguishable from the isotropic result while the energy loss with 9^ = is 
larger than the isotropic result for p^6 GeV. For example, for a 20 GeV bottom quark we find 
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FIG. 9: Anisotropic energy loss for = 1.0, 
Os = 0.3 and 9^ = (dotted line) compared 
to = 7r/2 (dashed line) function of 

velocity (normalized to the isotropic result, 
full line). 
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FIG. 10: Anisotropic energy loss for a bot- 
tom quark for ^ = 1.0, Og = 0.3 and 0^=0 
(dotted line), 9^ = tt/2 (dashed line) com- 
pared to the isotropic result (full line). 
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FIG. 11: Anisotropic energy loss for ^ = 
10.0, tts = 0.3 and 6^ = (dotted line) com- 
pared to 9 = tt/2 (dashed line) function 
of velocity (normalized to the isotropic re- 
sult, full line). 
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FIG. 12: Anisotropic energy loss for a bot- 
tom quark for ^ = 10.0, as = 0.3 and = 
(dotted line), 9^ = 7r/2 (dashed line) com- 
pared to the isotropic result (full line). 



that the energy loss at 9^ = is approximately 10% higher than the isotropic result. 

It should be noted that for very small velocities, the energy loss for the angle 9^ = becomes 
negative because of the presence of modes on the unphysical sheet as discussed in Sec. lIIII However, 
for ^ = 1.0 this occurs at very small velocities where the heavy quark approximation used is invalid 
in any case so that within the valid range of velocities {v > 0.39 and v > 0.71 for bottom and 
charm quarks, respectively) the energy loss is always positive. 

In Figs. ^2 and El we present the same plots but for ^ = 10. As can be seen from Fig. ^2 
for small velocities the energy loss is larger at 9jj = Tr/2 while for large velocities it is larger at 
9y = similar to the situation at = 1. From Fig. we see that the bottom quark energy loss 
at 9y = 7r/2 is less than the isotropic result while the energy loss at = is larger than the 
isotropic result for p^8 GeV. For example, for a 20 GeV bottom quark we find that the energy 
loss for 9^ = 7r/2 is approximately 10% lower than the isotropic result and for 9^ = the energy 
loss is approximately 20% higher than the isotropic result. Again for velocities less than v ~ 0.4 
the energy loss at 9^ = becomes negative; however, now the velocity where it becomes negative 
is on the border of the applicability of the heavy quark approximation for a bottom quark given in 
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FIG. 13: Anisotropic energy loss for ^ oo, 
Os = 0.3 and 9^ = (dotted line) compared 
to 9y = tt/2 (dashed line) as a function of 
velocity (normalized to the isotropic result, 
full line). 



FIG. 14: Anisotropic energy loss for a bot- 
tom quark for ^ — > oo, as = 0.3 and 6^ = 
(dotted line), 9^ = tt/2 (dashed line) com- 
pared to the isotropic result (full line). 



the previous paragraph. One is therefore lead to wonder what happens if we continue to increase 
the anisotropy. 

D. Infinite anisotropy limit 

As Figs. [TUl and [T^ demonstrate, as the anisotropy strength ^ is increased the directional depen- 
dence of the heavy quark energy loss increases. Additionally we found that the velocity at which 
the energy loss becomes negative increases as ^ is increased. It would be very instructive to then 
consider the limit of ^ ^ oo to see the limiting behavior of the energy loss. 

In Figs. El and El we present the same plots but for ^ = oo. As can be seen from Fig. 1131 for 
small velocities the energy loss is larger at 9^ = tt/2 while for large velocities it is larger at ^t, = 
similar to what we found at finite ^. From Fig. El we see that the energy loss of a bottom quark 
with 9y = tt/2 is less than the isotropic result while the energy loss of a bottom quark with 6y = is 
larger than the isotropic result for p^lO GeV; however, the results are not dramatically different 
than those obtained at ^ = 10. However, the velocity where the energy loss becomes negative, 
V ~ 0.65, is well within the range of applicability of the heavy quark approximation for a bottom 
quark (v > 0.39). We are therefore led to the conclusion that for very strong anisotropics heavy 
quarks with small velocities will experience an energy gain rather than an energy loss. Additionally 
we see that the energy loss is more negative at 9v = than any other angle of propagation. This 
means particles moving along the anisotropy direction experience the most energy gain. This is 
a rather surprising result but its origin can be traced to the presence of modes on the unphysical 
sheet as discussed in Sec. II I II These modes therefore represent a qualitative change to the physics 
in highly anisotropic plasmas. 

V. CONCLUSIONS 

In this paper we have calculated the complete leading-order collisional energy loss of a heavy 
quark propagating through an anisotropic quark-gluon plasma. The results were presented as a 
function of the heavy quark velocity and propagation angle. We presented a discussion concerning 
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the kinematical range over which the results obtained here apply when the quark mass is large 
but not infinite and we then explicitly applied the results to the energy loss of a bottom quark in 
an anisotropic setting. It was shown that for anisotropic systems the collisional energy loss has 
an angular dependence which increases as the coupling and/or anisotropy is increased as expected 
on intuitive grounds. Quantitatively, for Ug = 0.3 and a 20 GeV bottom quark we found that the 
deviations from the isotropic result were on the order of 10% for ^ = 1 and of the order of 20% for 
^ > 10. When translated into the difference between longitudinal and transverse energy loss this 
results in a 10% difference at ^ = 1, a 30% difference at ^ = 10, and a 50% difference at ^ = oo. 

More importantly, however, we found that for small velocities the sign of the energy loss becomes 
negative representing energy gain instead of loss whenever ^ > 0. In the isotropic limit and 
assuming an infinitely heavy quark it can be shown that the energy loss is positive definite and 
our results confirm this expectation for isotropic systems in contrast to the earlier calculation of 
Braaten and Thoma. In the anisotropic case, however, there is no guarantee that the energy loss 
will be a positive quantity. We have shown here that the negative contributions to the energy 
loss come from singularities on the neighboring unphysical sheets which are not present in the 
isotropic case For small anisotropics the velocities at which the energy loss becomes negative 
are lower than the thermal bound set by ~ \J'iT jra. For large anisotropics the velocities at which 
the energy loss becomes negative for a bottom quark are above the thermal lower bound and are 
therefore within the range of applicability of the infinite quark mass limit considered here. 

The results we have presented for the heavy quark energy loss should be compared with the 
expected momentum-space anisotropy generated during the early stages of RHIC and LHC heavy- 
ion collisions. At the earliest times that a particle description is appropriate the relevant scales are 
Vt ~ Qs where Qg is of order 1 GeV at RHIC and 2-3 GeV at LHC energies and ~ 1/r where 
T is the time at which the hard gluon occupation number drops below 1. In order to estimate 
T we follow the logic used in the bottom-up thermalization scenario giving parametrically 
QgT ~ as We therefore estimate that parametrically ~ aJ'^Q^ and as a result ^ ~ 
We therefore see that in the weak-coupling limit large anisotropics can be generated if the bottom- 
up scenario is true. Using a realistic coupling of a<j ~ 0.3 gives ^ ~ 6 but this estimate could 
change dramatically depending on the specific coefficients omitted in our parametric estimates. 
Additionally the time estimate Q^r ~ as ^^"^ is based on an implicit assumption of isotropy and 
collisional broadening and must be revisited in the context of an anisotropic quark-gluon plasma. 
Also in terms of producing phenomenologically relevant predictions it should be mentioned that 
since there are unstable modes present in the system, whatever value of one starts with, it is 
expected to go to zero very rapidly. In order to predict the total effect on observables we would 
need to fold together the results obtained here with the time evolution of ^ which is not currently 
available. 

It should be mentioned that thermal effects are only expected to make the energy loss more neg- 
ative at small velocities; however, it would be nice to have an explicit calculation of the anisotropic 
energy loss in the limit u — > with the infinite-mass assumption relaxed. Additionally, a re- 
finement of the calculation involving different momentum-space anisotropics for the quark and 
gluon distribution functions was suggested. Another area for improvement is that this work ig- 
nores next-to-leading order radiative energy loss. In isotropic systems it has been shown that the 
radiative energy loss of a heavy quark is larger than the collisional energy loss and thus cannot 
be ignored. As illustrated by this work, however, a calculation of the radiative energy loss for 
an anisotropic system would be considerably more involved than the equivalent calculation in an 
isotropic one. Nevertheless, it seems necessary in order to truly understand heavy quark energy 
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loss in an anisotropic plasma. 
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